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Abstract 

We present a one-parameter extension of the raise and peel one-dimensional growth 
model. The model is defined in the configuration space of Dyck (RSOS) paths. Tiles 
from a rarefied gas hit the interface and change its shape. The adsorption rates are 
local but the desorption rates are non-local, they depend not only on the cluster hit 
by the tile but also on the total number of peaks (local maxima) belonging to all 
the clusters of the configuration. The domain of the parameter is determined by the 
condition that the rates are non- negative. In the finite-size scaling limit, the model is 
conformal invariant in the whole open domain. The parameter appears in the sound 
velocity only. At the boundary of the domain, the stationary state is an adsorbing 
state and conformal invariance is lost. The model allows to check the universality of 
nonlocal observables in the raise and peel model. An example is given. 

1 Introduction 

In the study of adsorption- desorption (deposition-evaporation) processes on a planar sur- 
face, the Edwards- Wilkinson [1] and the Kardar-Parisi-Zhang [2] growth models have been 
extensively researched. The corresponding dynamical critical exponents z being equal to 

2 respectively 3/2. Recently the raise and peel model (RPM) was introduced [3]. In this 
one parameter dependent model, adsorption is local but desorption is not (the interface is 
"peeled"). When the parameter is changed in the critical domain, z changes too, varying 
from z = 1 to zero. The model was extended to contain sources at the boundaries [I] and 
defects [5]. The RPM was also used to check various estimators in information theory [6]. 
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If one fixes the value of the parameter such that z — 1, the model has magical properties. 
The system has a space-time symmetry (conformal invariance) and the stationary state 
probability distribution function (PDF) has fascinating combinatorial properties [7] being 
related to a two-dimensional equilibrium system [8]. This is also the single case where the 
system is integrable. 

In the present paper we will consider the RPM at the conformal invariant point only. 
With a few exceptions, it is not yet known which characteristics of the stationary state are 
universal and the time-dependent properties of the interface are by and large, uncharted 
territory. The aim of this paper is to introduce a new model, dependent on a parameter 
p, which is in the same universality class as the RPM. This implies that in the finite-size 
scaling limit the models coincide for the whole range of values of p. It turns out that in this 
limit, the single difference between the models is the value of the sound velocity v s (p) which 
is a function of p. Since the models are not local, finding models belonging to the same 
universality class is not an obvious matter. 

For reasons which will be apparent later on, we call the new model the peak adjusted 
raise and peel model (PARPM), it is presented in Section 2. Except for p — 1, the rates 
are dependent on the system-size L and on p. The PARPM for p ^ 1 is not integrable and 
therefore all the results we have are based on Monte Carlo simulations on large lattices. In 
the stationary states, the PDF's have no magic properties. 

In Section 3 we show that, in the finite-size scaling limit the density of contact points is 
independent of the parameter p. This observation is important since for p = 1 one recovers 
the RPM and for this case, the density of contact points is known exactly (this knowledge 
is based on conformal field theory and combinatorics [1]). 

We next examine the time dependence of the average number of clusters for various 
lattices sizes using the Family- Vicsek analysis [9]. The large time behavior of the density 
is given by a few levels of the spectrum of the evolution operator (Hamiltonian) and it is 
well understood for p = 1. We have observed that by a rescaling of the time, which implies 
changing the sound velocity v s (p) as a function of p, one obtains the same finite-size scaling 
function of tv s (p)/Lv s (l) for all values of p (L is the system size). In order to check that 
this picture is correct, we have done a finite-size scaling analysis of the spectrum of the 
Hamiltonians and found the same values for v s (p) as the ones determined in obtaining the 
Family- Vicsek function. 

The sound velocity stays finite for p in the interval < p < 2 and vanishes for p = 2 (the 
procedure how to take the limit p — > 2 is explained in the text). At p = 2 the stationary 
state is an adsorbing state and conformal invariance is lost. 

We believe that the results mentioned above are a clear demonstration that the PARPM 
is in the same universality class as the RPM. These observations opens the possibility to 
look for other universal quantities than those already considered in the study of the RPM. 
Since the models are not local, there is no recipe to find them. We will show that the average 
density of sites where desorption does not take place is universal: for large values of L, its 
value is independent of p. 

Our conclusions are presented in Section 4. 
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Figure 1: An example of a Dyck path for L = 14. There are four contact points and three 
clusters. The substrate profile is shown in blue. 

2 Description of the peak adjusted raise and peel 
model 

We consider an open one-dimensional system with L + 1 sites (L even). A Dyck path 
(restricted solid-on-solid (RSOS) configuration) is defined as follows. We attach to each site 
i non-negative integer heights hi which obey RSOS rules: 

h i+1 - hi = ±1, h = h L = (z = 0,l,...,L-l). (2.1) 

There are 

Z(L) = L!/(L/2)!(L/2 + l)! (2.2) 

configurations of this kind. If hj = at the site j one has a contact point. Between two 
consecutive contact points one has a cluster. There are four contact points and three clusters 
in Fig. [TJ 

A Dyck path is seen as an interface separating a film of tilted tiles deposited on a substrate 
from a rarefied gas of tiles (see Fig. [2]). The stochastic processes in discrete time has two 
steps: 

a) Sequential updating. With a probability P(i) a tile hits the site i — 1, . . . , L — 1 

(Ei-P(i) = !)■ In the RPM, P(i) is chosen uniform: P(i) = P = 1/(L-1). In the PARPM, 
this is not longer the case. For a given configuration (there are Z(L) of them), all the peaks 
are hit with the same probability R p = p/(L — 1) (p is a non-negative parameter), all the 
other sites are hit with the same probability Q c = q c / (L — 1) where 

q c = {L-l-pn c )/(L-l-n c ), c = 1, 2, . . . , Z{L). (2.3) 

Here n c is the number of peaks in the configuration (labelled by c). q c depends on the 
configuration c (with L sites and n c peaks) and on the parameter p. Obviously: 

n c R p + (L — 1 - n c )Q c = 1. (2.4) 

b) Effects of a hit by a tile. The consequence of the hit on a configuration is the same 
as in the RPM at the conformal invariant point. Depending of the slope = (/ij+i — frj_i)/2 
at the site i, the following processes can occur: 

1) Si = and hi > hi_\ (tile b in Fig [2]). The tile hits a peak and is reflected. 
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Figure 2: Example of a configuration with 4 peaks of the PARPM for L = 18. Depending 
on the position where the tilted tiles reach the interface, several distinct processes occur (see 
the text). 

2) Si = and hi < h^\ (tile c in Fig. [5]) . The tile hits a local minimum and is adsorbed 
{hi ^hi + 2). 

3) Si = 1 (tile a in fig. [2]). The tile is reflected after triggering the desorption (hj — > hj — 2) 
of a layer of b — 1 tiles from the segment {j = i + 1, . . . , i + b — 1} where foj > hi = h i+b . 

4) Sj = — 1 (tile in Fig. [2]). The tile is reflected after triggering the desorption (hj — > hj — 2) 
of a layer of 6 — 1 tiles belonging to the segment {j = i — 6+1, . . . , i— 1} where hj > hi = hi-b- 

The continuous time evolution of a system composed by the states a = 1,2,..., Z(L) 
with probabilities P a (t) is given by a master equation that can be interpreted as an imaginary 
time Schrodinger equation: 

j t Pa(t) = -J2H a , b P b (t), (2.5) 

where the Hamiltonian H is an Z(L) x ZiV) intensity matrix: H a b non positive (a 7^ b) and 
J2 a Ha,b = 0. —H ab is the rate for the transition — > \a). The ground-state wavefunction 
of the system |0), H\0) = 0, gives the probabilities in the stationary state: 

\0) = Y / P a \a), P a = \im P a (t). (2.6) 

a 

In order to go from the discrete time description of the stochastic model to the continuous 
time limit, we take At = 1/(L — 1) and 

Hac = -r ac qc (c^a), (2.7) 

where r ac are the rates of the RPM and q c is given by Eq. ( 12. 3p . The probabilities R p 
don't enter in (12 .5p since in the RPM when a tile hits a peak, the tile is reflected and 
the configuration stays unchanged. Notice that through the g c 's the matrix elements of 
the Hamiltonian depend on the size of the system and the numbers of peaks n c of the 
configurations. 

As an example, we consider a system with L = 6. In this case there are 5 configurations 
shown in Fig. [3J The Hamiltonian is given by: 
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Figure 3: The five Dyck path configurations for L = 6 
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recovers the Hamiltonian of the RPM for L = 6 





The unnormalized PDF in the stationary state is given by the eigenvector corresponding 
to the eigenvalue zero of (5): 



11(5 -P) 15(5 -p) 15(5 -p) 3(5 - p) \ 

2(5 - 3p) ' 4(5 - 2p) ' 4(5 - 2p) ' 5 - 2p ' / ' 1 ' J 

Except for p = 1 [3], there are no nice combinatorial properties of the components in 
f)2.9p for rational values of p. This was checked also for larger lattices. 

Let us observe that if p = 5/3, the configuration |1) which corresponds to the substrate, 
becomes an adsorptive state. This phenomenon is general. For an arbitrary value of L, the 
substrate has n p = L/2 peaks therefore, using (12. 3p we conclude that for p = 2{L — l)/L, 
the substrate is an adsorbing state. For larger values of p, one obtains negative rates, which 
gives 

<p < 2{L- l)/L (2.10) 

as the domain of p. 



3 Conformal invariance in the PARPM 

We have invented the PARPM expecting for p ^ 1 different properties than those observed. 
It came as a surprise that the PARPM belonged to the same universality class as the RPM. 
Let us first show that in the stationary state the density of contact points g(x, L) (x is the 
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Figure 4: The density of contact points g(x, L) divided by ( 13. ID for p = 0.01 and lattices 
sizes L = 510, 1024, 2048 and 4096. 



distance to the origin and L the size of the system) is, in the finite-size scaling limit (x » 1, 
L » 1 but x/L fixed), independent of p and equal to 

/ L \- 1/3 
g(x, L) = C(- sm(7rx/L)\ , (3.1) 

where 

x/3 

C = -—^T(-l/6) = 0.753149... . (3.2) 

67T6 

This expression is exact for p = 1 |4] and its functional form (not the constant C which 
is obtained using combinatorics) is a result of conformal invariance. In Figs. H] and [5] we 
show the x/L dependence of the densities of contact points divided by the expression (13. ip . 
as obtained from Monte Carlo simulations for two extreme values of p (0.01 and 1.99) and 
various lattices sizes. Similar calculations, done for other values of p, give the same result: in 
the finite-size limit, the density of contact points is independent of p. From the data collapse 
seen in the figures, we conclude that the first test of universality is successful. 

Let us consider now space-time properties of the model. One starts at t — with the 
configuration describing the substrate and we look at the average number of clusters k(t, L) 
as a function of time and lattice sizes. In order to determine the dynamical critical exponent 
z for various values of p, we compute the Family- Vicsek scaling function [9] 

K(t,L) = k(t,L)/k(L) -1, (3.3) 

where k(L) is the average number of clusters in the stationary state. In the finite-size scaling 
limit one expects 

K(t, L) = K{t/L z ). (3.4) 



(3 




Figure 6: The Family- Vicsek scaling function K(t,L) as a function of t/L for p = 
0.1,0.5, 1, 1.5 and 1.9 and lattice sizes L = 1024,2048,4096 and 8192. 
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Figure 7: The Family- Vicsek scaling function K(tv s (p), L) as a function of tv s (p)/v s (l)L for 
p = 0.1, 0.5, 1.5, 1.9 and L = 1024, 2048, 4096, 8192. The 20 sets of data collapse on a single 
curve. 

In Fig. |6]we show the results of the Monte Carlo simulations. For each value of p one has 
data collapse which implies z = 1 for all the values of p (see (13.4p ). One now show that by 
a change of the time scale, one can get the different scaling functions in Fig. [6] to coincide. 
Let us keep in mind that in a conformal invariant theory, the finite-size scaling limit of the 
spectrum of the Hamiltonian has the following behavior (the ground-state eigenvalue E is 
equal to zero for a stochastic process): 

lim EAL) = irv s Xi/L, i = 0, 1, 2, . . . (3.5) 

L— >oo 

where v s is the sound velocity and Xi are critical exponents. This implies that 

K(t/L) =J2°i eM-Eit) =J2°i exp(-7iv sXi t/L), (3.6) 

i i 

where the constants Cj depend on the initial conditions. As discussed in [10] , for p — 1 one 
has: 

3 

v s (l) = -V3, x x = 2, x 2 = 3, ... . (3.7) 

For large values of t/L one has obtained a very good fit to the data using only x\ and x 2 - A 
p-dependent change of the time scale would correspond to a change of the sound velocity. 

In Fig. [7] we show the data of Fig. [6] choosing for v s (p)/v s (l) the values 1.538 (p = 0.1), 
1.300 (p = 0.5), 0.703 (p = 1.5) and 0.472 (p = 1.9). Notice that the 20 sets of data collapse 
on a single curve. 

In order to check that the determination of the sound velocity is correct, we have di- 
agonalized the Hamiltonians up to L = 26 sites for p = 0.5, 1 and 1.5 and looked at the 
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second energy gap E 2 (p). The calculation of E 2 (p) is easier to calculate numerically as 
compared with E\ since it corresponds to the smallest eigenvalue in the parity odd sector 
of the Hamiltonian. We have computed the ratios £ , 2(0.5)/i?2(l) and E2{1.5) / E 2 {V) which 
should converge for large L to v s (0.5)/v s (l) respectively v s (1.5)/v s (l). Extrapolants give 
the values 1.307 respectively 0.6972 in excellent agreement with the values obtained from 
Family- Vicsek scaling. 

We believe that in this way, we have shown that the PARPM belongs to the same 
universality class as the RPM for all values of p. 

We give now an example which shows why it is useful to have models in the same 
universality class. Let us consider, in the stationary states, the average density of sites in 
which the slop vanishes (local minima and maxima): 

r(L) = l/(L-l)X;(l-| ai |). (3.8) 

i=i 

This quantity (ADSSV) was extensively used in the study of the phase diagram of the RPM 
|10j . There is a conjecture [IT] about the values of r(L) at the conformal invariant point of 
the RPM: 

t{L) = (3L 2 -2L + 2)/(L - 1)(4L + 2). (3.9) 

This conjecture was checked on various lattice sizes. For large values of L, r(L) has the 
following behavior: 

t(L) = A + B/L + 0(1/L 2 ), (3.10) 

where A = 3/4 and B = —1/8. Are these coefficients universal? Using Monte Carlo 
simulations for lattices of different sizes (for L up to 65000) we obtained the following results: 
p = 0.5, A = 0.7500,5 = 0.219; p=l,A = 0.7500,5 = -0.125; p = 1.99, A = 0.7500,5 = 
-0.302. 

We conclude that the average density of sites where the slope vanishes is a universal 
constant equal to 0.75 with the leading correction term being non-universal. An explanation 
of this interesting result in the context of conformal invariance is missing. It is also not clear 
if the ADSSV is a universal observable in other models defined on Dyck paths. 

One can use the asymptotic value of t(L) to have an estimate of the sound velocity. The 
average value of the density of peaks, equal to half the value of r(L), is < n c > av = 3/8. 
Substituting this value in (12 .4p . and using ( 12 .7p one obtains 

< q c (p) >av= —ttt = — =— . 3.11 
v s (l) 5 

We have assumed that the fluctuations of the density of sites where the slop vanishes are 
small. One can compare the values obtained for v s (p) using (12 .4p with those obtained from 
the diagonalization of the Hamiltonian and from Family- Vicsek scaling. In Fig. [8]we compare 
some of the results obtained from the Family- Vicsek scaling with the prediction in (13. lip . 
The agreement is excellent. One notices that if p approaches the value p — 2, v s (p)/v s (l) 
approaches the value 2/5. For the substrate the density of peaks is equal to 1/2, this gives, 
using (12. 7p the value v s (2) = for p = 2. As a consequence, v s (p) decreases from the value 
8/5 for p = to the value 2/5 Hp approaches the value 2 from below and has a discontinuity 
at p = 2. At p = 2 the substrate is an absorbing state. 
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Figure 8: Comparison of the ratios v s (p)/v s (l) given by (13. lip (red line) and some of the 
estimated values (black dots) from the Family- Vicsek finite-size scaling. 

4 Conclusions 

The raise and peel model was up to now the single known conformal invariant growth model. 
The model described in this paper is different: it has rates which are dependent on global 
aspects of the configurations (the number of peaks). The fact that conformal invariance is 
maintained came as a surprise. The new model is not integrable and the stationary state 
does not have the remarkable properties of the RPM which are a consequence of integrability. 
Conformal invariance was tested for several values of the parameter p by computing the 
density of contact points, the spectrum of the evolution operator and the Family- Vicsek 
scaling function. In the finite-size scaling limit, only the sound velocity depends on p. 

Using the new model we have discovered that the average density of sites where the slope 
vanishes is, for large lattice sizes, a universal quantity and equal to 3/4 for any value of p. It 
would be interesting to look for other quantities which are universal. The model described 
here can be used tool to check universality. 

As explained in the text, the parameter p varies in the domain < p < 2(L — 1)/L. For 
each value of p < 2, one considers values of L > 2/(2 — p) which makes sure that the rates 
are positive. We have observed that if p approaches the value 2, the sound velocity decreases 
to a limiting value v s = 2/5, and one observes a slowing down in the time dependence of 
physical processes as one approaches the stationary state (see Figure 6). 

An interesting phenomenon occurs if one chooses p = 2(L — 1)/L. Conformal invariance 
is lost and the stationary state is an absorbing state which is the substrate (see Figure. 1)). 
Although the stationary state is unique, for lattices sizes L > 90, the relaxation time is very 
long and increases exponentially with the size of the system. One observes metastable states. 
If the initial condition are changed, the system ends up in different metastable states. A 
detailed presentation of the p = 2(L — 1)/L case will be published elsewhere [T2] . 
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